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Abstract 

In this paper random walks on the Penrose lattice are investigated. 
Heat kernel estimates and the invariance principle are shown. 



1 Introduction 

The Penrose tiling [9l [15] is the most famous nonperiodic tiling of the plane 
and is an unfailing source of beautiful properties and phenomena to be ex- 
plored [1]. There is associated with the Penrose tiling by the usual duality 
relation a lattice that may be called the Penrose lattice. Kunz [HJ discussed 
conditions under which simple nearest-neighbor random walk on the Pen- 
rose lattice would conform to a version of the central limit theorem, but 
did not completely resolve the issue. In the present paper, it is proved 
that random walk on the Penrose lattice satisfies the invariance principle 
(in annealed sense, see detailed explanation in the next section) that with 
appropriate scalings, the random walk process converges weakly to a non- 
degenerate rotation-invariant Brownian motion. Szasz [17] has conjectured 
that a related Lorentz scatter system also satisfies this principle. The proof 
of the discrete, random walk counterpart of the conjecture, also formulated 
by Szasz, is the main result of the paper. 

Central to the approach followed here is the notion of roughly isometric 
weighted graphs, defined precisely later, which share diffusion properties (cf. 
[To]). This enables us to relate the simple random walk on the Penrose 
lattice to translation invariant walks on the plane square lattice 1? . Let 
y) denote the graphical distance (number of edges in the shortest path) 
between vertices x and y. Many random walk processes obey two-sided 
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Gaussian heat kernel estimate {GEa^2) 

exp ( j <Pn{x,y) < ^^exp I ) (1^ 



for < n < d{x, y), with a > 1, c > and C > 0. Here pn = Pn+Pn+i, where 
Pn{x,y) is the probabihty that a walker departing from vertex x is found at 
vertex y after n steps. Delmotte |3] has shown that the bounds {GE2): 

c [ Cd{x,y f\ \ ^ ( cd{x,yy 

exp < Pn [x, y) < — ^ exp 



V (x, ^/n) \ n J ' V" (x, y/n) \ n 

(2) 

(where V (x, r) denotes the volume of geodesic balls in graph distance) are 
stable under rough isometry for random walks on weighted graphs (which we 
define precisely later) . In other words if two graphs are roughly isometric and 
([2]) holds for one then it holds for the other as well. We know that ^ holds 
for the simple symmetric random walk on ( V" (x, r) ~ r^, i.e. satisfy ([T]) 
with a = 2) and consequently holds for the Penrose graph if it that graph is 
roughly isometric to Z^. A very short, direct proof will be given of rough 
isometry between the Penrose graph and 1?. 

Unfortunately the exponents in the upper and lower estimates in (P and [2]) 
contain different constants which reflects local inhomogeneities. The rough 
isometry invariance and the Gaussian estimate ([2]) have been proved along 
a series of estimates in which some cumulation of constants is unavoidable. 
Consequently if we are looking for the central limit theorem or for the invari- 
ance principle we need a different approach. 

In the field of stochastic processes and statistical physics a powerful 
method has been developed to investigate random walks in random envi- 
ronment, on percolation clusters and interacting particle systems. A key 
result in this direction is the celebrated paper by Kipnis and Varadhan [12] 
and its influential extension by De Masi, Ferrari, Goldstein, Wick [H] in which 
annealed central limit theorem and invariance principle is shown, (the initial 
environment is averaged with respect to an invariant measure), (for further 
details see Section H] and [16] on convergence notions in random environ- 
ments). That result provides us immediate derivation of the central limit 
theorem and the invariance principle for the random walk on the Penrose 
lattice in the same sense. 

In what follows we introduce the basic terminology then the statement is 
proved. 
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2 Preliminaries 



We will consider infinite connected graphs with vertex set T. If an edge joins 
vertices x and y we denote that edge by x ~ ?/. The distance d (x, y) will 
be the shortest path metric. In particular we will speak about the integer 
lattice Z*" = (Z™, d) graph where vertexes are elements of Z*" and x,y & Z™ 
form an edge, x ~ ?/, if and only if |x — ?/| = 1. We will speak about the 
integer lattice (Z™, |.|) if we consider the same vertex and edge set but the 
metric is the Euclidean one. We do not define the Penrose tiling ( cf. [5]), 
we assume that it is well defined and given for us on M^. 

The Penrose lattice (P, |.|) is a metric space. It is the set of centers 
(centroids) of the tiles equipped with the Euclidean distance. Two tiles are 
neighbors if they are edge adjacent. Two vertexes of the Penrose lattice are 
neighbors if they centers of neighboring tiles. Those vertexes form edges of 
the lattice. We will speak about Penrose graph, with the same vertex and 
edge set but with d (x, y), the shortest path graphs distance. Let P = (P, dp) 
denote the Penrose graph, and {l?,di) integer lattice graph. 

We distinguish tilings by fixing a reference vertex and identifying it with 
the origin of M^. We denote by f2 the space of tilings (union of ten tori 
fij.j C i, j = 0..4, i < j , for details cf. Section 2.1 in [11]) • Let d (x) be 
the degree of x, the number of neighbors. 

Definition 1 A graph is weighted if a symmetric weight function fi^y > 
is given on the edges. This weight defines a measure on vertexes and sets: 

(1 (x) = ^ ] fJ-x^y 

y^x 

lx{A) = 5^/i(x) 

We denote the ball of radius of r by B (x, r) = {y : d (x, y) < r} and we call 
V (x, r) = fi{B (x, r)) its volume. In particular for the Penrose lattice (and 
graph) fi^^y = 1 if x y is an edge, while n^^y is zero otherwise. The same 
applies for the integer lattice. 

A Markov chain, with transition probabilities P (x, y) is reversible (/z- 
reversible) if there is a /i measure such that ji (x) P (x, y) = fi (y) P {y, x). 

Definition 2 In general a random walk Xn on T, a weighted graph, with /x 
is a reversible Markov chain defined by the one step transition probabilities: 

P{X„, = |/|X„_i =x) = P (x, y) = 
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The random walk on the Penrose lattice (and graph) is reversible Markov 
chain with transition probability P{x,y) = l/d{x) = 1/4 for x ~ y. It is 
clear that d (x) P (x, y) = d (y) P {y,x) = 1. Denote Xj the actual position of 
the Markov chain (random walk) which is well-defined for any fixed Xq G F 
and it is the reference vertex of the Penrose lattice. 

3 Heat kernel estimate for the Penrose graph 

First of all we give the definition the bi-Lipschitz property and rough isom- 
etry. 

Definition 3 A metric space (F, d) is bi-Lipschitz to (F', d') if there is a 
bijection $ from F to T' and a constant C > 1 such that for all x ^ y eT 

^d {x, y) < d' ($ (x) , $ {y)) < Cd (x, y) (3) 

Definition 4 Two weighted graphs F with fi and F' with fi' are roughly iso- 
metric (or quasi isometric) (cf. /S", Definition 5.9]) if there is a map (j) from 
F to V such that there are a,b,c, M > for which 

-d (x, y)-b<d'{4> (x) , (y)) < ad (x, y) + b (4) 
a 

for all x,y eT, 

d'{<p{V),y')<M (5) 

for all y' G F' and 

-/i (x) < jj! ((f) (x)) < cfi (x) (6) 
c 

for all X G F. 

Remark 1 It is clear that if from F to T' is a rough isometry then there 
is a rough isometry 0' from F' to F as well. . 

Theorem 1 (Solomon [14^ ) The Penrose lattice is bi-Lipschitz to the integer 
lattice. 

Proposition 2 The Penrose lattice is rough isometric to the integer lattice. 
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The statement follows from the bi-Lipschitz property. A very short and 
direct proof can be given, which we present here. 

Proof. [Proof of Proposition [2] Denote by m the smaller distance between 
the opposite boundaries of the thin rhombus and write e = A/2m/4. Consider 
the integer lattice eZ^. It is clear that if an open square with edge length 
e contains a center of a rhombus that it is fully contained by the closed 
rhombus. Let \Ef map the center of the rhombus to the center of the square. 
It is clear that is rough isometry from the Penrose lattice to the integer 
lattice . ■ 

Proposition 3 The Penrose graph is roughly isometric to the integer lattice 
graph. 

Proof. Let us consider \I', the map introduced above, between F and eZ^. 
Now we consider the graph distances dp, d^. It is clear that 

dp {x,y) < d^{^{x),^{y)). 

The opposite inequality is also easy. Let 2L be the maximal number of 
squares which is needed to cover the largest diagonal of rhombi. It is clear 
that L is bounded since the diameter of the rhombi is also bounded. Then 

c?z(^(x),^(y)) <2Lc/p ix,y). 

It is also clear that the conditions (ISlTG]) are satisfied. ■ 

Lemma 4 The Penrose lattice and the Penrose graph are roughly isometric. 

Proof. Let $i the rough isometry from the graph F to the lattice Z^, let $2 
the identity map on which is bi-Lipschitz between the integer lattice and 
graph, finally $3 the rough isometry from the integer graph to the Penrose 
graph. The existence of $3 follows from Proposition |3] and Remark [1] Then 
$ = $3 o $2 o $1 is rough isometry between the Penrose lattice and graph. 
■ 

Now we recall Delmotte's result [4] omitting the third equivalent state- 
ment, the parabolic Harnack inequality, since we do not need it in the sequel. 

Theorem 5 Let F with fi be a weighted graph. Assume that there is a po > 
such that for all edges P {x, y) > po. Then the following statements are equiv- 
alent. 
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1. there are C, c > 0, a > 1 such that for all x,y and n > 

I r^ d{x,yf \ ^ C ( d{x,yy 

■ exp -C < pn [x, y) < — ^ exp -c- 



V (x, y/n) \ n I ' y {x, y/n) \ n 

(7) 

holds, 
2. 

{%) The volume doubling condition iyD) holds: there is a C > such that for 
all X & r , r > 1 

V{x,2r) < CV (x,r) 

and 

(a) the Poincare inequality {PI2) holds: 

there is a C > 0, such that for all x and, r > l,f : i? (x, r) — )■ M 

E (/(?/)- /b)V (y) < E (/(i/)-/(^))V,,. 

y(iB{x,r) y,z(iB{x,r) 

where fB= v^)^yeB f iv) iv) ^ / ^ 0' B = B{x,r). 

It is well known that the volume doubling property as well as the Poincare 
inequality are rough isometry invariant. Of course volume doubling can be 
replaced with V (x, r) ~ r" and ([7]) reduces into ([T]) . 

Corollary 6 If T with /i and T' with ^' are roughly isometric graphs then 
{GEa,2) holds for one if and only if holds for the other. 

Theorem 7 The Gaussian estimate (Gi?2,2) holds for the random walk on 
the Penrose graph. 

Proof. Proposition [2] ensures that Penrose graph is rough isometric to Z^. 
It is well-known that (G'-E'2,2) holds for the random walk on the integer lattice 
(graph), and then by Corollary [6] ((^-^2,2) holds for the random walk on the 
Penrose graph as well. ■ 



4 The invariance principle 

In this section we confine ourself to the Penrose lattice. The Gaussian esti- 
mate (G£^2,2) provides a nice description of the random walk on the Penrose 
graph but the different constants in the exponents mean that we have only 
estimate of the variance and it may change from place to place as well as 
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in time. Particularly we do not know if the properly scaled mean square 
displacement {(P (Xo,X„)) has a limit. Of course we expect that due to 
the asymptotic spherical symmetry of the Penrose tiling the diffusion matrix 
is the identity matrix up to a fixed constant multiplier. In other words the 
scaled mean square displacement is direction independent. In order to obtain 
the invariance principle for the Penrose lattice we need a different method. 
This is the method of ergodic processes of the environment "seen from the 
tagged particle" (cf. [12], [6]). Thanks to the result of De Masi & all [6] it is 
enough to check that the conditions of Theorem 2.1 in [6] are satisfied and 
that the covariance matrix is positive definite. 

There are several formulation of the invariance principle, (see for the 
classical formulation in We say that Xt satisfies the central limit theorem 
(CLT) if there is a cr > such that 

in distribution. What is slightly stronger, the process can be re-scaled, that 
is for £ ^ for all t 

X(^) = ^ (t) (8) 

where Wiener process with variance at in the sense of finite dimensional 
distributions. The invariance principle holds if the convergence holds for the 
path-space measures for the processes (cf. Goldstein [S]). This requirement 
is equivalent with ([8]) and "tightness" of the process (cf. [2]). 

Let the path-space measure of the processes (xi^^ ) ,6>0. Follow- 

\ / t>o 

ing [6] we say that X*^^^ = (x^'^^) converges weakly in u-measure to Y if 

V / t>o 

for any continuous function F on the path-space D ([0, oo) , F): 

{F (x(^))) ^E^(r). 

Definition 5 We say that the invariance principle holds if the weak conver- 
gence in ^-measure to the Wiener process holds in this sense. 

We consider the environment process a;„ seen from the particle. It is 
more convenient to use (as it is done by Kunz in [TT]) the Markov chain 
Zn = {un,Xn). Let US uotc that (F (X^^))) means that the underlying 
environment process is started from the invariant measure /i, E^ (F (X*^^))) = 
E^ (F (X*^^^) \ujq = lo) oj is chosen according to /i, in other words we average 
with respect to the (initial) environment and obtain an annealed type of 
result (cf. [16]). 
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Theorem 8 The random walk on the Penrose lattice satisfies the invariance 
principle with non- degenerate covariance matrix. 

Proof. Kunz has shown that Zn is ergodic (see also a more general result 
by Robinson [13j). The invariant measure is combination of the Lebesgue 
measure A on the tori Qij. For < z < j < 4 

where r is the " Golden mean" (a/5 + l) /2. It is clear that X„ = Yli=i ^ (^i-i? ^-i 
where V Zi) = Xj— is an antisymmetric function, (cf. [H] (2. 3), (2. 6) 
and the remark below it.) It follows from the definition that X„ and Zn as 
well are reversible. The random walk is well defined for all u, for all Pen- 
rose lattice, with given reference vertex, hence we have the path metric P^j- 
Similarly is well defined if the initial lattice is chosen according to the 
invariant measure. Let denote the corresponding expected value. The 
only properties are to check that the conditional drift 

</. = E^(Xi-Xo|Xo) (9) 

exists and that the covariance matrix 

D = E^((Xi-v^)(Xi-<^)*) (10) 

is non-degenerate. For any given Xq = x the conditional drift evidently 
exists thanks to the bounded distances of neighbors. The conditions of the 
main result of [B] are satisfied, hence the invariance principle holds for Xi in 
the sense of Definition O 

Let us recall that the D always exists (see Remark 1. below (2.30) in 
0). 

We show that the covariance matrix is positive definite. Let us con- 
sider the annulus B (0, C2y/n) \B (0, Cia/ti) intersected with the cone about 
a given direction e G with angle a fixed 7r/2 > a > 0. Let H denote 
the intersection. The constants Ci, C2 are arbitrary and fixed. Let us recall 
(Lemma Hj) that the Penrose lattice and graph are roughly isometric 

^-E {e*X^X:Ee\X,=x,) 



n 

^E((e*X„)'|X„ =xo) 



n 



> (exf Pn{xo,x) 



n 



\H\ 2 c" 

> c' (cos (a) \/n) — exp 

n n 



^ (aC2v^ + &) 
n 



> 
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independently of Xq and e, hence the covariance matrix is non-degenerate. 
(Here the effect of all previous constants are absorbed into the last constant 
c.) By this we have shown that the invariance principle holds for the random 
walk on almost all Penrose lattice and the limiting process is a non-degenerate 
Brownian motion. ■ 

Remark 2 The results presented in this paper carry over easily to other 
quasicrystals which can he constructed by the projection methods similar to 
the one produces the Penrose tiling. This applies to generalized Penrose 
tilings (produced by p- grids), higher dimensional Penrose tilings and stochas- 
tic tilings. 

Remark 3 It seems plausible that with some extra work one can show that 
the covariance matrix is the identity matrix multiplied with a positive con- 
stant. The exact value of the constant ought to be determined as well. 
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